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RENEWAL SEQUENCES, DISORDERED POTENTIALS, 
AND PINNING PHENOMENA 

GIAMBATTISTA GIACOMIN 

^—>. . Abstract. We give an overview of the state of the art of the analysis of disordered 

^s^ ' models of pinning on a defect line. This class of models includes a number of well known 

and much studied systems (like polymer pinning on a defect line, wetting of interfaces 
on a disordered substrate and the Poland- Scheraga model of DNA denaturation). A 
remarkable aspect is that, in absence of disorder, all the models in this class are exactly 
r~^, , solvable and they display a localization-delocalization transition that one understands 

^Sj ' in full detail. Moreover the behavior of such systems near criticality is controlled by 

a parameter and one observes, by tuning the parameter, the full spectrum of critical 
behaviors, ranging from first order to infinite order transitions. This is therefore an ideal 
set-up in which to address the question of the effect of disorder on the phase transition, 
notably on critical properties. We will review recent results that show that the physical 
,S^ • prediction that goes under the name of Harris criterion is indeed fully correct for pinning 

T^ ' models. Beyond summarizing the results, we will sketch most of the arguments of proof. 
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1. Pinning and disorder: models and motivations 

1.1. The basic example: pinning of simple random walks. It is somewhat custom- 
ary to introduce pinning models by talking of pinning of simple random walks (SRW). 
This is due to a number of reasons, like the widespread grasp on SRW, or the fact that 
modeling several pinning phenomena naturally leads to random walk pinning, as we will 
see. However, we will see also that, in a sense, the SRW case is the hardest to deal with: 
nonetheless, we are going to follow the tradition and start from SRW pinning. 

Let S := {«S'„}„=o,i,... be a sequence of random variables such that 5*0 = and such 
that {Sn — Sn-i}n=i,2,... are IID (i.e., independent and identically distributed) symmetric 
random variables taking only the values +1 and —1. The disorder is given by a sequence 
cv := {uJn}n=i,2,... of real numbers and we will play with two real parameters /? and h. 
We actually assume that a; is a realization of an IID sequence of standard Gaussian 
variables (see Remark 11.51 for some comments on generalizations) . We call P the law of 
a; and we denote by 9 the left-shift operator on M^: {9uj)n = iOn+i- Our aim is to study 
the probability measure Pn,uj {N is a positive integer: we will be interested in the limit 
N ^ oo) defined as 
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N 



PN,ujiso, si,...,sn) ■■= ^^^— exp V] {(3uJn + h) l^„=o P(so, Si, ... , sn), (1.1) 

^^'- V^ J 

where 

(1) P{si,S2, ■ ■ ■ ,sn) = (1/2)^ if and only if sq = and \sn — Sn-i| = 1 for n = 
1,2,. ..,iV; 

(2) Zj\f^^ is the normalization constant (partition function), that is 

Zn,uj = ^ exp y^ (PuJn + h) lg„=o P{so, Si,..., sn)- (1-2) 

so,si,...,Sff \n=l / 

We will actually prefer a slightly different definition of the model, namely given the 
sequence s = {so, si, . . .} we set 

^^(s) := ^ exp (y1 il^^n + h) l.„=o J . (1.3) 

'"^ \n=l / 

Notice that this time Pn,uj is a measure on (infinite) sequences, namely the trajectory 
of the walk all the way to infinity, while in (jl.ip we had defined a measure only up 
to step (or time) N. As a matter of fact, if we consider cylindrical events of the type 
E = {s = {s„}„=o,i,... : sq = to, si = ti, . . . , Sn = tN}, then the measure of E under Pn,oj 
defined in (|1.3p coincides with PN,ojito, ti, ■ ■ ■ , tN)- 

Remark 1.1. It is worth stressing that, unless /3 = 0, in this model there are two sources 
of randomness: the polymer chain is modeled by a random walk with law P and the 
disorder is a typical realization of the random sequence w with law P. These two sources 
of randomness are treated in very different ways: u is quenched, that is chosen once and 
for all, while the polymer location fiuctuates and in fact we study the distribution of S 
under Pnuj- 
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Figure 1. Two trajectories sampled from Pn,ui for different values of /3 and h, with A*' 
very large, and represented as directed polymers, in the sense that we plot the function 
n H^ Sn or, equivalently, we look at the walk {(n, S'„)}„=o,i,... with one deterministic 
component. The case A sketches a delocalized trajectory that is what one observes when 
h is, for example, negative and large (one should first think of the homogeneous case 
/3 = 0). Since the SRW is a periodic Markov chain, the origin is visited only at even 
times, so only u>2,^4,^6, ■ ■ ■ play a role and they are the only charges (a charge at n is the 
quantity /3aj„ + h) marked in the drawing, with a filled circle if they are visited and with 
an empty circle if they are not visited; not all the unvisited charges are marked). The 
distinctive aspect of case A is that there are just a few contacts, i.e. visits to the origin, 
and then the walk resembles a walk conditioned not to hit zero. The case B is instead 
what one observes when h is positive and large: the number of contacts is large, as a 
matter of fact the drawing wants to suggest that there is a positive density of contacts. 
It is natural to call such a regime localized, in contrast to the previous one that we call 
delocalized. Note that both case A and B are atypical for the free walk {P = h — 0) in 
which there is a zero density of contacts but they are spread through the system and the 
walk certainly does not stay on one side of the axis as in case A. It is important to remark 
that we have a full up-down symmetry and this implies that in the delocalized regime A 
the walk is either delocalized above or below the axis with probability 1/2. Of course at 
this stage it is highly unclear that one observes either localized or delocalized trajectories 
(for typical ui) and to a certain extent this is not correct because one has to exclude the 
so called critical regime, which however appears only at exceptional values of l3 and h 
(the phase transition point, or critical point). We have of course avoided the delicate 
issue of the role of disorder, at the hearth of this presentation. Here we will simply 
content ourselves with pointing out that, for example, even when h is very large and 
negative (pushing thus toward delocalization) any amount of disorder, i.e. /3 > 0, yields 
a positive density of sites in which l3ij„ + h is positive and therefore attractive. There 
could therefore be a smart targeting strategy of the polymer in placing the contacts at 
these sites, leading thus to localization. 



As it is well known, the Markov process S is null-recurrent, namely every site of the state 
space Z is visited (infinitely often) P-almost surely, but the expectation of the time between 
successive visits to a given site is infinite. Let us be more explicit about this last concept 
and let us introduce, for m G Z, the random variable Ti(m) := inf{n > : Sn = m} and. 



for j > 1, also Tj{m) := inf{n > Tj 



-i[m 



Sn = ITT,}. It is then a direct consequence of 
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the (strong) Markov property that {Tj^i(m) — Tj{m)}j=i2^... is a sequence of IID random 
variables. It should be also clear that the law of {Tj+i{m) — Tj{m)}j=i^2,... does not 
depend on the value of m: we are going to denote Tj(0) simply by tj and, since Sq = 0, 
we are setting tq = 0. The recurrent character of S boils simply down to the fact that 
^^ P(ri = n) = 1 and the recurrence is of null type because E[ti] = +oo: the distribution 
of Ti is known in detail [23^ Ch. Ill] and in particular 

P(n=2n)"^°^-^, (1.4) 

where c = l/yivr and we have introduced the notation a„ ~ 6„ for lim„_too cLn/bn = 1- 
Since (clearly) E[ri] = oo, the classical Renewal Theorem (see (II. lip below) tells us that 
the expected number of visits to of S up to time N is o{N) (a more precise analysis 
shows that it is of the order of viV , see |23l Ch. Ill] or Theorem 11.21 below). 

As we shall see, the trajectories of the process S are very strongly affected if (3 or h 
are non zero and, except for critical cases, what happens is roughly that, in the limit as 
N ^ oo, under Pat^o; the expected number of the visits paid by 5" to is of the order of 
N, or it is much smaller than y/N (in some cases one can show that they are 0(1)). A 
first glimpse at these different scenarios can be found in Figure [H 

1.2. The general model: rene'wal pinning. We have introduced the r sequence in the 
previous subsection in order to give some intuition about the model, but its interest goes 
well beyond. A look at (jl.Sp suffices to realize that Z^v^^ can be expressed simply in terms 
of r: 



^N,LU 



Eexp ( ^(/3u;„ + /i)l„g^ J , 



(1.5) 



where we have introduced a notation that comes from looking at r = {tj}j=o,i,... as a 
random subset of N U {0}, so that n £ t means that there exists j such that tj = n. 
Therefore the model in (II. 3p is just a particular case of when r is a general discrete 
renewal: 



rlP* 



r) 



J — exp ^ {(iu}n + h) Iner , (1-6) 

'V:'-' \n=l / 



dP ^ ^ • Zj, 

where the superscript f, that stands for free, has been introduced because a slightly 
different version of the model is going to be relevant too: 



dP^,. 



(r) : = 



1 / ^ \ 

^— exp ^ {(3L0n + h) InGT 'i-N^r, (1-7) 

^,<^ \n=l J 



dP ^ ' z% 

and c stands for constrained. Let us stress that by (discrete) renewal process r we simply 
mean a sequence of random variables with (positive and integer valued) IID increments: 
we call these increments inter-arrival variables. 

We will not be interested in the most general discrete renewal, but we will rather focus 
on the case in which 

K(n) :=P(Ti=n) = ^, n G N = {1, 2, . . .} (1.8) 

where a is a positive number and 

lim L{n) = c^ > 0. (1.9) 

n— >oo 
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We call K{-) inter-arrival distribution. Note that we always assume K(0) = 0. Moreover 

we will assume that that Y2neN ^i''^) — ^'- *^^ ^^^^ YlneN^i''^) "^ -'- ^^^ *o interpreted 
as the case of a terminating renewal, in the sense that K{oo) := 1 — ^„gj^ K{n) is the 
probability that ri = +cxd. Therefore, if K{oo) > 0, the cardinality \t\ of the random set 
r is almost surely finite. The case of K{oo) = is instead the case of a persistent renewal, 
and \t\ = oo almost surely. But persistent renewals are of two different kinds: they are 
positive persistent if X]n'^-^("')(~ Eti) < cxd, or null persistent if the same quantity 
diverges. This terminology reflects the fact that for any aperiodic renewal (aperiodicity 
refers to the fact that ri does not concentrate on a sublattice of N) the law of large numbers 
ensures that almost surely 

lim -|rn[0,n]| = -^ G [0,1], (1.10) 

n— >oo n Eti 

so that if Eri = cxd we are facing a zero density renewal. We stress that the last statement 
holds also for terminating renewals for which E[ti] = '^n(zfqnK{n) + oo-R'(oo) = oo (note 
on the way that, for us, ^„ . . . never includes n = oo). 

Very relevant for the analysis of renewal processes is the renewal function n i-^ P(n G r), 
that is the probability that the site n is visited by the renewal. We call such a function 
K{-)-renewal function when we want to be more precise. The asymptotic behavior of the 
renewal function is captured by the so called Renewal Theorem (for a proof see e.g. [8]). 
This theorem says that if r is an aperiodic renewal (the generalization to the periodic case 
is immediate) we have 

lim P(neT) = -^ G [0,1]. (1.11) 

n^oo Eti 

Note the link with (jl.lOp . but note also that this is little informative if Eri = oo. The 
leading asymptotic behavior in such a case is summed up in the following statement that 
calls for the deflnition of the Gamma function: T{x) = J^ t^~^ exp(— t) di, for x > 0. 
Recall that in our set-up E[ti] = oo either because i^(oo) > (terminating renewal), 
regardless of the value of a, or because a < 1. 



Proposition 1.2. Assuming (jl.Sp and (|1.9p we have: 
(1) ifK{oo) >0 then 



P(neT)"~°"^^. (1.12) 



K{ooy 



(2) // K{oo) = andae (0, 1) then 



n-*oo /asin(7ra)\ ^_-^ 



Proposition 11.2( 1) is a classical result detailed for example in [9| or [29l, § A. 6]. Propo- 
sition [L2l!|2) is instead more delicate and while the case a S (1/2, 1) is under control since 
|27j . the full case is instead a rather recent result [21j . Note however that if full proofs 
are non-trivial, one can rather easily flnd intuitive arguments suggesting the validity of 
Proposition O [29., § A. 6]. 

The asymptotic behavior for the case a = 1 is known too, but this case is a bit anomalous 
and, in this review, we will often skip the results for a = 1 that would make the exposition 
heavier. 
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Remark 1.3. In full generality, given a renewal r one can find a Markov process 5 on a 
state space containing a point such that ri = inf{n E N : S'n = 0}. Still in full generality, 
the state space can be chosen equal to N U {0}, see [23 App. A. 5] for details. Therefore, 
with this remark in mind, one could go back to the original definition ()1.3p without loss 
of generality. 

Remark 1.4. Everything we are going to present works assuming simply that K{-) is 
regularly varying or, equivalently, that L(-) is slowly varying. Examples of slowly varying 
functions include log(n)'^, c a real number, or any product of powers of iterated logarithms 
(see [9] for full definitions and properties or [29^ § A. 4] for a quick sum-up). Regularly 
varying functions are a natural set-up for pinning models also because some natural cases 
do involve slowly varying functions: for example the law of ri for the two dimensional 
simple random walk one has K(n) ~ c/(n(logn)^), for an explicit value of c > |40j . 

In this subsection we have focused on the behavior of the free system: /3 = /i = 0. The 
observations made in the caption of Figure [1] do apply to the general case too (that is to 
{P,h) 7^ (0,0)), with, nevertheless, two distinctions: 

(1) if a > 1 and if the renewal is persistent, the free renewal is already localized, since 
by the Renewal Theorem the contact points have a positive density. We will see 
that this affects the discussion in the caption of Figure [1] only for what concerns 
the critical case, but the general picture still holds; 

(2) if r is terminating it is of course harder to localize the process, but in reality it is 
rather easy to show that the model can be mapped to a persistent r case, precisely 
if one sets K{n) = K{n)/{1 — K{oo)) and if r is the (persistent!) iC(-)-renewal 

Z^,^ = E exp(^(/3cu„ + /i + log(l-if(oo)))l„e^J ;iVG? , (1.14) 

and the same is true also at the level of the measure Pn,uj itself. The proof of such 
a statement is absolutely elementary and it is detailed for example in [291 Ch. 1]: 
note that, from a mathematical standpoint, this allows us to restrict ourselves to 
persistent renewals r. 

1.3. A gallery of applications. The localization mechanism captured by class of models 
we have just introduced comes up in modeling a variety of phenomena. Here we just extract 
some examples and cite some references. 

Polymers and defect lines. The interaction between polymers, chains of elementary units 
called monomers, and the surrounding medium or other polymers is omnipresent in physics, 
chemistry and biology. We cite for example [28] , but it is of course impossible to account 
for the literature in such a direction. The case we are interested in is the one in which a 
polymer is fluctuating in a neutral medium except for a line, or a tube, with which the 
polymer interacts. Actually, also cases in which the line is for example a surface or even 
simply a point may be modeled by the type of pinning models we are considering. The key 
point is that polymers are often modeled by self-avoiding random walks and a simplified 
way to impose the self-avoiding condition is considering directed walks (see Figured]). So 
the polymer pinning model becomes precisely the renewal pinning we are considering (we 
refer to |17l [Ml [261 B2l [57] for examples of the phenomena that are modeled via directed 
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walk pinning). Here we just stress that the dimensionahty of the problem enters the re- 
newal pinning only via the exponent a: for example a polymer in three dimensions pinned 
to a line can be modeled by {(n, 5'n)}n=o,i,...) where 5" is a random walk in two dimensions, 
for which a = (see Remark 1 1.4p . The general case of a physical space of (i+ 1 dimensions 
leads to a = ((i/2) — 1, for d > 2, and of course a = 1/2 if d = 1. 

Wetting 'phenomena. Modeling interfaces in two dimensional media by random walks has 
a long history [1] that is somewhat summed up also in [29] App. C] in which one can 
find the explanation of why anisotropic Ising models do reduce in a suitable limit to the 
renewal pinning model with q = 1/2. A particular choice of the boundary conditions 
leads to the so called wetting problem |131 l46j. which is just the case in which the random 
walk trajectories that one considers are only the ones above (and touching) the axis: with 
reference to Figure [H to obtain allowed trajectories one has to flip over the negative 
excursions. As it is explained in detail in |29t Ch. 1], this problem just corresponds to 
renewal pinning with a = 1/2 and K{po) = 1/2 and, at the level of contact points, the 
process can be mapped to the case in Figure [U with h replaced hy h — log 2 (see (I1.14p ). 

DNA denaturation: the Poland-Scheraga model. Two-stranded DNA has been often mod- 
eled by two directed walks with pinning potentials, see e.g. |47j and references therein. 
Since the difference of two independent random walks is still a random walk, we are back 
to renewal pinning. However directed walk models lead to values of a that are in contrast 
with observation, so that a considerable amount of work has been put into understanding 
whether (in our language) renewal pinning is a reasonable model and which a should be 
chosen (see in particular [51], but once again we refer to [29] for a more complete bibli- 
ography). We note that inhomogeneous or disordered modeling is really required in this 
context, because the pinning strength does depend on the type of base pair, see Figure HI 
The renewal pinning model with inhomogeneous charges has been and is extensively used 
for the study of DNA denaturation [11^117): appropriate values of a are close to 1.15. 



Loops 




T6 - T5 



Tg -rs 



TU - Tl3 



Figure 2. The two thick hnes are the DNA strands. They may be paired, gaining thus 
energetic contributions that depend on whether the base pair is A-T or G-C (the model 
is therefore inhomogeneous). There are then sections of unpaired bases (the loops) to 
which an entropy is associated. The DNA portion in the drawing corresponds to the 
renewal model trajectory with tj — tj-i = 1 except three r-interarrivals (so the loops 
correspond to inter-arrival of length 2 or more). 



Remark 1.5. For DNA denaturation taking uji Gaussian is not appropriate. In this case 
uJi should rather be a binary variable. To be more precise one should also take into account 
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stacking energies, that is energies depend on blocks of two pairs, and probably one should 
also study correlated sequences of bases. Sticking to the issue of binary variables versus 
Gaussian ones, we take this occasion to stress that much of the mathematical literature is 
written for rather general charge distribution (say, with finite exponential moments of all 
orders). 

2. The homogeneous case 

The full solution of the non disordered (/? = 0), or homogeneous, case is crucial and, at 
the same time, it is rather elementary once it is phrased in the renewal theory language. 
Such a solution has been repeatedly presented in the physical literature in several par- 
ticular instances, by using what a probabilist would call generating function techniques: 
in particular one can find a very nice and complete presentation in [23]. However the 
presentation we are going to outline in detail here is different and much more direct. 

In this section, but also later, with abuse of notation we will denote by Z'^ j^ (a = c, f ) 
the partition function Z'^ ^ when (3 = 0. Let us start by observing that we can write 

N n 

Zkh = J2 Y. U'^W^Kii,). (2.1) 

n=l ^eN": j=l 

Note that if /i = then Z^ ^.^ = P{N £ r), i.e. the partition function is just the K{-)- 
renewal function. The right-hand side of (j2.ip is still a renewal function if e^K{-) is an 
inter-arrival law. And indeed it is if X^„gi^ e K{n) < 1 and in this case Z^ ^ is the e K{-)- 
renewal function: its asymptotic behavior is hence given in Theorem II. 2 1, but we prefer to 
delay such a result since a unified approach holds for every h. In fact if X]„gj^ e^K{n) > 1 
we can renormalize the expression by introducing an exponential correction, going back to 
a renewal function (times an exponentially growing factor). Precisely we call 6(> 0) the 
(unique) solution of 

^exp{-bn + h)K{n) = 1, (2.2) 

neN 

and we set Ki,{n) := exp(— 6n -|- h)K{n). We have therefore defined a function h h^ b{h) 
for h such that X^neN ^'^-^C") — ^' ^^^* ^^ ^•^'^ ^ — ^c(O), with 

/i,(0) := -\ogY,K{n). (2.3) 

n 

For h < hc{0) we set instead b{h) = and Ko{n) := exp{h)K{n) (of course the latter 
notation is poor since h is not explicit). With this notations we can write 

N n 

Z^^;, = exp(6iV) J] Yl llKhil,) = expibN)Ph{NeT), (2.4) 

n=l £gN": j=l 

where, under P/j, r is a Kuh\{-)-Tenewal. By the Renewal Theorem 

hm P,(iV G r) = -i-^, (2.5) 

which is a positive constant if /i > /ic(0), but it is zero if /i < /ic(0) because the Kb(/i)(-)- 
renewal is terminating. For h = /ic(0) this limit may or may not be zero, but let us 
postpone this issue to Remark 12.11 Let us focus for now on the fact that for h < /ic(0) the 
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Renewal Theorem does not yield the leading behavior, but thanks to Propositioi jl.2l fl) we 
see that 

N^oo K{N) 

(l-exp(/i)(l-K(oo)))' 
With these results in our hands we see 



PhiN G r) "~°° ,,,:^ ^ ^^, . , (2.6) 



(1) that since A^ ^ logP/j(A^ G r) vanishes as A^ — > oo, we have therefore proven that 

F{0,h) := Jim ^logZ^,;, = b{h), (2.7) 

for every h. The quantity f(0, /i) is usually called free energy (per unit volume) 
and of course we use such a notation because later there will be f(/3, h); 

(2) that (j2.4p goes well beyond Laplace asymptotics: this is very relevant and allows 
us for example to compute the limit of 



^N,h (ti = ni,T2=ni+n2,...,Tj =ni-\ h Uj) = J|(e''i^(nj)) 



^N-ni-...-nj,h 



yc ' 

(2.8) 
as A'^ ^ oo. For example when h > /ic(0) the ratio of partition functions in the 
right-hand side converges to exp(— (ni + . . . + nj)F(0,^)) and therefore the all 
expression converges to 11^=1 -^F{o,/i)('^i)- It is rather easy to see that the same 
holds also for h < /ic(0). 

Remark 2.1. In the above list we have been a bit clumsy about the critical case h = hc{0), 
but in reality what happens at h = hc{0) is clear too. First of all, X^„-ft^F{o,hc(o))('^) = 
1, so that the associated renewal is persistent. More precisely ^f(o,/ic(o))(") = ^i') if 
^„ir(n) = 1, and Zf^ ^ = P{N £ r), and otherwise -ft^F{o,/ic(o))(") is just a multiple of 
K{-) and Z^ ^ coincides with the Kp(o,/i^(o))(0"'^6newal function computed in A^. Recall 
now that the KprQfitQ-j\{-)-renewal function converges to a positive constant if a > 1 and 
to zero otherwise. But when it converges to zero there is Proposition II. 2l f2) that comes to 
our help so that once again we know the sharp asymptotic behavior of Z^^. In particular 
F(0,/ie(0))=0. 



All these remarks are telling us in particular that (j2.2p is a formula for the free energy, 
in the sense that f(0, /i) = 6 if there exists a positive solution b to (|2.2p and otherwise 
f(0, /i) = 0. From such a formula one can extract a number of consequences that are 
summed up in the caption of Figure [3j In particular the behavior of the free energy near 
criticality is trivial for h < /ic(0), but it is not for h > /ic(0): let us make explicit the 
behavior of f(0,/i) as /i \ /ic(0). If X]n'^^(^) ^ °° ^^^ ^^ hc{0) = (which we may 
assume without loss of generality: recall p.l4p !l 

^{1 - eM-Hh)n))K{n) = 1 - exp{-h) '"^^ h. (2.9) 

n 

The asymptotic behavior of the left-hand side is easily obtained since for every fixed n 
the limit of (1 — exp(—b{h)n)/b as 5 \ is n. On the other hand 1 — exp(— x) < x for 
every x > 0, so that the Dominated Convergence Theorem yields that the left-most side 
in (12. 9p is asymptotically equivalent to b^^nK{n), and therefore b{h) ~ h/^^nK{n). 



10 
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f(0,/i) 




dhmh) 




Figure 3. The function h h^ f(0, /i) is non decreasing, convex and non negative (con- 
vexity follows either from (|2.2|l or (|2.7|l '). It is therefore equal to up to hc{0) — sup{h : 
f(0, h) = 0} and after this point it is positive and strictly increasing. Of course h — hc{0) 
is a point of non analiticity of the map h ^^ F{0,h): this map is of course analytic on 
(—oo,hc{0)). It is analytic also on (/ic(0),oo), by the Implicit Function Theorem for 
analytic functions. The graph of dhP{0,h) indicates that we are considering the case 
Q = 1/2: f(0, •) is C^ but not C^. For a G (1/2, 1] the slope of /i h^ dhF{0, h) at KiQ) 
is infinite and for a > 1 a jump discontinuity appears. We stress that dh'P{Q, h) is the 
contact fraction of the system, see Remark 12.51 and therefore such an observable has a 
jump at the transition for a > 1. 



If instead a G (0, 1) formula ()2.9p still holds, but the asymptotic behavior of the left-hand 
side is gotten by Riemann sum approximation: 



^ (1 - exp(-.(ft)„)) A-(„) "^» ,'C,,^ ^-^j-f" 



(bn) 
h-c.l i^^S^dx = 6"c.i:(i^, (2.10) 



nl+a 



a 



h\0 



SO that it suffices to invert the asymptotic relation (&(/i))°Cx(r(l — a)/a) ~ h and this 
of course gives that b{h) is asymptotically proportional to h^''^. 

The arguments that we have developed directly lead to the following statement (see 
Ch. 2] for a more complete statement and for more details on the proof): 



Theorem 2.2. The critical behavior of the free energy is given by 

(5\o I ci6 if a > I, 



F(0,/ic(0) + <5) 



with 



ci := 



l-K(oo) 



and C2 




(2.11) 



(2.12) 



E„eN"^(^) ' V Cii-r(l-a) 

Moreover in full generality, as N ^ oo, P^/j (that denotes the measure P^^^ when (3 = 0) 

converges weakly in the product topology ofM.^ to a limit measure P/j. The limit process 
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is a KpfQi^\{-) -renewal, namely: 

k 
Ph{n=£i,T2 = ei+£2,...,rk = ei + ...+£k) = llK^io,h){£j), (2.13) 

i=i 

where 

jeMh-nFiO,h))K{n), z/f(0,/i)>0, 
'^^(°''^)^"^ = \eMh)Kin), ^fHO,h)=0. ^''''^ 

Therefore f(0, h) > implies that the limit process is positive persistent (in fact, the inter- 
arrival distribution decays exponentially), while instead if f{0, h) = the limit inter-arrival 
distribution has power law decay and, if h < /ic(0), the K^iqi^\[-) -renewal is terminating. 

Remark 2.3. Once sharp results for the constrained case are obtained, one can deduce 
sharp results on the free case. This is just based on the elementary formula 

N 

Z^,, = Y,K,hK{N-n), (2.15) 

n=0 

that we have written in the case X^„gN K{n) = 1 and we have introduced 

K{N) := Y^ K{n). (2.16) 

nGN:n>Af 

For example \i h> hc{0) from (|2.4p we have 

N 

Zj^^^ = exp(F(0, h)N) Y, Ph{n G r) exp(-F(0, h){N - n))K{N - n). (2.17) 

n=0 

Since exp(-F(0, /i)(iV - n))K{N - n) is bounded below by Yl,j>N-m^^v{-'P{'^^h)i)K{i) 

and since X^„=o-^('^ ^ t)K{N — n) = 1 for any persistent i^(-)-renewal we obtain that 
for every A'' 

Zj^j^ > exp(-/i)exp(F(0,/i)iV). (2.18) 

A (rough) bound in the other direction is obtained by neglecting P/i(n G t)K{N — n) in 
the right-hand side of (I2.17p . so that 

Z%,h < i_p(o;,) ^^P(F(O.MA^)> (2.19) 

which holds once again for every A^. The sharp asymptotic result is 

„f N^oo exp(F(0, /l)Af) ^^^ / /n i.\ \-77t \ 

^"'^ ~ E,.»A-„..,(») g;''^''^'°-''"'^'"' 

^ (l-exp(-/.))a.P(0/.) 
1 — exp(— f(0, nj) 

This type of estimates directly leads to computing the limit behavior of P^/j, see [291 
Ch. 2] for details. 
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Remark 2.4. A key concept in statistical mechanics (and a key concept here) is the 
notion of correlation length. For example a natural correlation length of the system for 
h > hc{0) is given by the rate of exponential decay, as n — > oo, of P/i(?^ G ''") to its limit 
value l/E/i[ri]. One can show [30J in particular that if h is sufficiently close to /ic(0) then 
Ph{n e r) - l/Eh[n] > and 

hm - log (Phin G r) - -^) = -f(0, h), (2.21) 

which says that the correlation length coincides with 1/f(0, h). Even if one takes a Enite 
volume viewpoint, 1/f(0, h) appears as a natural correlation length, for example because 
in ()2.4I) one sees that is only when N is of the order of 1/f(0, h) that one starts observing 
the exponential growth of the partition function. One could push these arguments a bit 
further and see that if A^ is much smaller than 1/f(0, h) (of course this has a precise sense 
only when h \ hc{0)) then P^ ^ resembles P, while for N much larger than 1/f(0, h) the 
measure P^ ^ starts exhibiting localization. The fact that the inverse of the free energy 
is the correlation length still holds also in presence of disorder |351 [521 [53] , even if a full 
understanding of this important issue is still elusive. 

Remark 2.5. The density of contacts, or contact fraction, that is the limit as A^ ^ co 
of A^~^E,j[X;^=i InGr] coincides by Theorem [22] with lim„^ooP^(n G r) = l/E^ri. Note 
moreover that N~^'Eh['^^^^ Ingr] = N~^dhlogZf^ i^, so that the contact fraction is equal 
to dhF{0,h) (except, possibly, at h = /ic(0)). Moreover, by simple conditioning arguments 
one easily sees for example that for h ^ /ic(0) 

1 



lim lim max 

J\/^oo N-*oo n:M<n<N-]\I 



Pk,i(" e r) 



and, by arguing like in Remark 1 2. 3 [ one directly sees that he same statement holds for the 
free case. 

3. The disordered case 

3.1. The quenched free energy. An elementary observation that turns out to be really 
crucial for us at several instances is that for every M = 0, 1, . . . , A^ 

log Z^,, > log Zl,^^ + log Z^_M,eM.- (3-1) 

It is simply proven by restricting the renewal trajectories, in the expression for Z^^, to 
the ones that contain the contact site M and by using the renewal property. By averaging 
over the disorder one sees that {ElogZ^^}7v=o,i,... is super-additive and this entails [l3] 
the existence of the limit 

hm ^ElogZ^_^ =: f(/?,/i), (3.2) 

and the limit of this sequence coincides with its supremum: 

F(/3,/i) = sup^ElogZ^^^. (3.3) 

NeN ^v 



The super-additive property (|3.ip can however be exploited further in order to obtain re- 
sults about the limit of the non averaged sequence {log ^^^}Ar=o,i,..., by applying the tools 
available for super-additive ergodic sequences, and notably the celebrated Kingman's The- 
orem [43] . Alternatively one can stick to the super-additive character of {E log Z^ ^j}n=q,i,... 
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and establish a concentration property of the non averaged sequence, either by using con- 
centration inequaUties, e.g. [45^ I51j. or even by more elementary tools [29' Ch. 4]. In all 
cases the result that one obtains is the existence and the self-averaging character of the 
free energy of pinning systems: 

Theorem 3.1. The sequence < N^^logZf^ ^> converges to f(/3,/i) both P(da;)- 

almost surely and in the L^ sense. 

Remark 3.2. It is not difficult |29t Ch. 4] to show that for every K(-), (3 and h there 
exists c > such that 

^^,. < ^^,. <ciVZ^,^, (3.4) 

uniformly in uj. We can therefore restate Theorem 13.11 replacing the superscript c with f . 

Another elementary central fact is that 



Zf^^ > E 



and therefore 




exp V(/3a;„ + /i)l„6, ; r n [1, iV] = {N} 



exp{l3ujN + h)K{N), (3.5) 



F{P,h) > 0. (3.6) 

We now partition the parameter space of the system into: 

C := {(/3, h) : f(/3, h) > 0} and V := {(/3, h) : f(/3, h) = 0} . (3.7) 

C and T> stand respectively for localized and Pelocalized regime, a nomenclature that 
calls for further explanations (see § 13.21 just below), but for the moment we just point out 
that one of our main aim is to characterize these regions as precisely as possible. And a 
substantial help is given by the fact that the function {(3, h) >-^ f(/3, h) is convex, as limit of 
convex functions, and it is monotonic non decreasing in both variables (monotonicity in h is 
immediate, in /3 it is instead a consequence of convexity and of the fact that 5/jE log Zf^ ^ = 
for /3 = 0). Since by (13. 6|) we see that V coincides with {(/3, h) : f(/3, h) = 0} so that V 
is a convex set. 

One can go beyond: Jensen inequality {annealing) yields 

exp((/i + /?V2)X]l„er);iVGr , 



n=l 



(3.8) 



ElogZ^,^ < logEZ^^^ = logE 

so that 

F(/3,/i) < F(0,/i + /3V2), (3.9) 

and if we recall that f{(3, h) > f(0, h) we directly get 

/ir"(/3) := hM - ^ < KiP) < h,{0). (3.10) 

As a matter of fact the upper bound can be made strict, that is /ic(/?) < ^c(O) as soon as 
/3 > 0, in full generality (the generality here refers to the choice of K{-), [6j) and in the 
framework that we consider here one can show also that, given K{-), for every /Jq > one 
can find an explicit constant c G (0, 1/2) such that /ic(/3) < /ic(0) - 0/3^ for /3 e (0, /3q] [291 
Ch. 5]. Instead, showing that hc{(3) > /i^""(/3) is a more delicate issue (and it is not true 
in general!). These bounds are summed up in Figure [Hand they imply that, since 2? is a 
convex set, then hc{-) is concave and, since it is bounded, it is continuous. 
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Jo)^-....,,^^,^^^,^^^^ y^ 




7"x>\ 7^'^^"^----,_ 









/ir"(/^) := ^c(O) - /3V2 / \^V \. 




hc{P) ^'\^ \. 



Figure 4. The critical curve (3 \~-* hc{l3) that separates D and £ is concave decreasing. 
This follows from the fact that I> is a convex set and from the explicit bounds we have 
on the critical curve. The upper bound h{P) is less explicit then the lower bound, 
but we stress that /i(/3) < hc{0) for every /3 > and this shows that disorder may 
induce localization and it never suppresses it. The lower bound comes from the standard 
annealing procedure. Note that the annealed partition function ¥.Z%^^, a = c,f , is just 
the homogeneous partition function with pinning potential h + /3^/2. 

Let US sum up the outcome of the arguments we have just outhned: 

Proposition 3.3. // we set hc{[3) = inf{/i : f(/3, h) > 0} then hc{f3) = sup{h : f(/3, h) = 
0} and 

C = {{p,h) : h>h^{f3)} and V = {{f3 , h) : h < hcif3)} (3.11) 

Moreover the function (3 i— > hc{(3) is concave, decreasing and (j3.10p holds for every (5. 



Remark 3.4. From (jS.Sp we actually extract the important observation that localization 
can be observed in finite volume, in the sense that (/?, h) ^ Cii and only if there exists N 
such that ElogZ^^ > 0. 



3.2. On path behavior. Characterizing localization and delocalization simply by looking 
at whether the free energy is positive or zero may look, from a mathematical standpoint, 
rather cheap. This is not the case as one can first see by observing that 

N 



dh^\ogZ%^^ 



E 



N,uj 



El 

.n=l 



n&T 



(3.12) 



which, by exploiting the convexity of the free energy and Theorem 13.11 tells us that P( do;)- 



a.s. 



dhFip,h) 



lim E 



N~ 



N,ui 



1\- 



n=l 



(3.13) 



when 9/jF(/3, h) exists. By convexity, such a derivative exists except possibly at a countable 
set of points and in any case (j3.13p can be extended to a (standard) suitable statement 
also if the derivative does not exist, in terms of right and left derivatives ^29j (as a matter 
of fact, in [35] it is shown that f(/?, •) is C°° except possibly at /ic(/3)). In the end (13.130 
is telling us that the contact fraction, i.e. the right-hand side in ()3.13p . of our system is 
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zero if /i < hdP) (that is, in the interior of 2?) and it is positive if h > hc{(3) (that is, in 
the whole of C). By itself this fully justifies our definition of (de)localization. 

However (|3.13p is still a poor result and plenty of questions could be asked about the 
limit of the sequence {P^^}Ar=o,i,...i a = c or a = f , starting with the existence of such a 
limit. And of course the question is: how close one can get to the very sharp description 
of the limit measure available for homogeneous systems? 

Work has been done in this direction, but we will not concentrate on this aspect. We 
just point out that 

(1) The localized phase is, to a certain extent, rather well understood. In the sense 
that if (/3, h) £ C than one can show that the weak limit as A'^ tends to infinity 
of the sequence of probability measures {Pn,u}}n exists P(da;)-a.s. and the limit 
process is a point process with a positive density of points [35]. One can show also 
other estimates going toward the completely clear picture that emerges from the 
homogeneous case. Intriguing differences however do arise, naturally connected 
to the existence of exceptional deviations in the sequence of charges. Moreover a 
number of open questions still stand (see e.g. [36J). 

(2) Progress has been made only recently on the delocalized phase, at least away of 
criticality [32] (see [52] for some estimates at criticality) . Essentially one now knows 
that in the delocalized non critical regime the number of contacts for a system of 
size N is 0(logA'^) and such a result has been achieved by a subtle argument 
combining concentration bounds and super-additivity properties of log Z^^. Such 
a result still leaves open intriguing questions in the direction, for example, of 
the precise results proven in [161 [39] in the homogeneous or weakly inhomogeneous 
context, see for example in the bibliographic complements at the end of [,29', Ch. 8]. 

3.3. The role of disorder. The main questions we want to address are: 

(1) How does the disorder affect the phase diagram? Namely can we determine hc{(3), 
for /3 > 0, beyond the bounds in Figure HP 

(2) What can one say about the critical behavior of the free energy? This amounts to 
estimating how f(/3, h) vanishes as h\ hdP). 

It is particularly interesting to raise such questions because we know /ic(0) and we know 
the sharp asymptotic behavior of f(0, h) for h close to hc{f3) (see Theorem 12. 2p . so that in 
our framework inquiring about the role of the disorder makes perfect sense. At this point 
it is important to underline that such questions do find partial (non rigorous) answers 
in the physical literature: the rest of this subsection is devoted to explaining what one 
expects on the basis of formal expansions, following some renormalization group ideas. We 
must say that the arguments that follow are adaptation to the pinning model context of an 
argument developed by A. B. Harris [37j in the context the Ising model with random bond 
defects. Harris' argument is based on the idea that the behavior of a system near criticality 
should become rather independent of fine details, so in particular one can replace the 
system by a coarse grained one without changing substantially the properties. What one 
actually tries to do is defining a renormalization transformation, like decimation or block 
summation, that, once applied repeatedly at criticality, transform the system into a limit 
model. Harris' work aims at determining whether introducing the disorder modifies the 
fixed point of the renormalization transformation: if the renormalization transformation 
suppresses the disorder and the limit point is like in the homogeneous case, then one says 
that disorder is irrelevant. If instead disorder is enhanced one says that disorder is relevant 
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and most probably the renormalization transformation flow leads to a fixed point which is 
different from the one obtained in the homogeneous case. It should be noted on one hand 
that at the border between relevance and irrelevance the renormalization transformation, 
to first order, neither decreases nor increases the disorder: this is the so called marginal 
case. On the other hand, Harris argument is just a small disorder expansion and as such it 
does not apply to the whole range of the parameters and, above all, it does not characterize 
the limit fixed point when disorder is relevant. 

Harris' ideas have been first applied in the pinning model context by G. Forgacs, J. M. 
Luck, Th. M. Nieuwenhuizen and H. Orland [25] and then by B. Derrida, V. Hakim and 
J. Vannimenus [20] with predictions that differ somewhat in a sense that we are going to 
explain just below. 

Let us start with an expansion that is freely inspired by [25] . Without loss of generality 
we assume /ic(0) = (recall (I1.14p ). Moreover the argument does not feel the boundary 
condition: we work it out in the free case. In what follows 6 := h + /3^/2 > 0: this change 
of variable is a natural one because in particular 



TV 



exp (5 ^ 1„6^ 

V n=l / 



E [ZijJ = E 
We set Cn = exp(/?a;n — /3^/2) — 1 and let us note that 
Elog 



zi 



N,5- 



(3.14) 



7* 



ElogEi^ 



ElogE^, 



N 



exp^(/3w„-/3V2)l„g^ 

\ra=l 

■ N 



i-ner) 



n=l 



= Elogl l + J];CnP^,5(^e^)+ Yl CniCn2P^,5({ni,n2}Cr) + ...j . 
\ n ni<n2 / 

(3.15) 

Let us now expand the logarithm and let us use the fact that the C random variables are 
centered and IID with variance equal to exp(/?^) — 1 to see that 

N 



Elog 



7^ 



i(exp(/32)-l)^P^,(nGr)2 + ... 



(3.16) 



n=l 



F(/3, /ir"(/5) +d) = F(/?, 6 - (3^2) = F(0, 6)-- (exp(/32) - i) (dsF{0, 6)f + . . . (3.17) 



By Remark 1 2. 5 1 for 5 > and as long as n and N — n are large, P^^(n G r) is close to 
dsF{d,5) so that from (|3.15|) we extract 

1 

2 

Of course this expansion is only formal and in order to make it rigorous one has to control 
the rest. Let us note on the way that one can in principle try to compute all the terms 
in this expansion, but the issue of controlling the rest is still there and convergence issues 
may very well require (3 to be small (note that we are expanding using as small parameter 
the variance of (", but aiming at capturing the critical behavior, hence h is small too). All 
the same, (I3.17P is compatible with hc{P) = hc{0) if f(0, (5) vanishes much slower than 
{dsF{0,6))'^ as (5 \ (/? possibly small, but fixed). But by Remark 13.131 (or directly by 
taking the h derivative in ()2.2p ) we see that dsF{0, 6) = l/E^ri and by direct computation 
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(similar to (j2.10p ) one sees that dsF{0,6) ~ (c2/a)5~^^^'° for a € (0, 1) (c2 is given in 
Theorem 12.21 but the precise value does not play a role here), while the contact fraction 
is bounded away for zero when a > 1 even approaching criticality. So (j3.17p is compatible 
with hciP) = hciO) if 

^l/a -^^O 0-1+1/a) ^^ ^ < 1. (3.18) 

This argument therefore suggests that disorder is irrelevant for a < 1/2. 

If a > 1/2 the expansion we have performed looks hopeless, but we may argue that this 
is just due to the fact that hc{P) > h'^'^{(3) and we are expanding around the wrong point. 
Of course we do know that f(/3, hc{(3)) = and therefore (|3.17p suggests that for /? small 
the shift of the quenched critical point 5c{l3) := /ic(/?) — h'^'^^{P) is found by equating the 
two terms in the rightmost side of (|3.17p and this procedure suggests Sc{P) ~ ^2o/(2a-i). 

A second approach is instead inspired by [20j. If we aim at analyzing whether the 
annealed system is close to the quenched system one could sit at the annealed critical 
point {h = hc{0) — /3^/2 = — /3^/2, i.e. 6 = 0) and study the variance of Z^ ^ (once again, 
the argument would go through also with constrained boundary condition). Divergence 
of the variance, as A^ ^ oo, would be a sign that quenched and annealed systems aren't 
close. Since at 5 = we have EZ^^ = 1 and 



varp(Z^^) = E 



(^U -1 



EE 



®2 



exp V(/?w„ - /3V2)(lner + Iner') - 1 



(3.19) 

with T and r' independent copies of the same renewal process. Integrating out the to 
variables we obtain 



N 

varp {Z%^^) = E' 



exp P^ ^ Inernr' - 1 



=1 



(3.20) 



This expression can be evaluated in a sharp way because the random set r n r' is still a 
renewal process and therefore the variance that we are evaluating is the partition function 
of a homogeneous pinning model (minus one). And the first relevant question is whether 
r n r' is a terminating or a persistent renewal. The inter-arrival law of r n r' can be 
expressed in terms of the inter-arrival law of r only in an implicit way, but the renewal 
function of r fl r' is explicit in terms of the renewal function of r: 

P®2(nGTnr') = P(nGr)^ (3.21) 

and T n r' is terminating (respectively, persistent) if ^^ P (n G r) < oo (respectively, 
Y^^ P (n G r) = oo) and, by Proposition 11.21 we see that 

oo ^ ^ 

72 :=5]P(nGr)2<oo ^ ^ _^_^ < oo ^ a < -. (3.22) 

n=l n 

By the general solution of the homogeneous model, cf. Section [21 we see that if r n r' is 
persistent, then for every /3 > the variance of Zj^ ^ grows exponentially, while if r n r' is 
terminating then X := ItPit'I — 1 is a geometric random variable (this is just a consequence 
of the renewal property) of expectation 72, that is P(A = n) = (72/(1 + 72))"(1/(1 + 72)), 
n = 0, 1, . . .. Therefore, as long as /? < /3o := •y/log((l + 72)/72), with p2 := (1/(1 + 72)) 



18 GIAMBATTISTA GIACOMIN 



we have 

Jim varp (Z^, J = , _ (, _^^') ^^^(^2) " ^ = 72/?^ + ■ • • (3.23) 

where the expansion is of /3 small. Therefore if rPlr' is terminating (note that r and r' are 
persistent since we are assuming /ic(0) = 0) the variance of Zf^ ^, at the critical annealed 
point, stays bounded and it is small if /3 is small. To complement (I3.23P note that 

supvarp Z^,^ < ^' -^^ - 1 < ^/J^, (3.24) 

N 1 - (1 -P2)exp(/3^) 

for some c > 0. 

Let us sum up the outcome of the arguments we have just outlined: 

(1) Both approaches suggest that disorder is irrelevant if a < 1/2 (and, as a con- 
sequence, relevant if a > 1/2, with the case a = 1/2 as marginal one), as one 
can read from (j3.18p and (|3.22p . Moreover the arguments do suggest that the 
annealed system is very close to the quenched one, in particular hc{(3) = hc{0), 
at least for /? not too large. This observation may be considered as the Harris 
criterion prediction for pinning models. 

(2) There is a difference between (I3.18P and (I3.22p in the case a = 1/2 that we cannot 
appreciate since we are assuming (II. 9p . In the more general framework of Re- 
mark [T3] one sees that the fact that L(-) diverges at infinity does not imply that 
T n r' is terminating, while it is sufficient to conclude that f(0,(^) is much larger 
that {dsF{0, (5))^ for 6 small. As a matter of fact, we are dealing with the marginal 
case in the renormalization group sense. This is a very subtle issue, still unresolved 
even on a purely heuristic level. We should stress that the steps that we have just 
presented here are just a part of the arguments in [25j and [20], in particular [25] 
aims at an expansion to all orders and [20] contains a subtle attempt to study the 
renormalization group flow for 5 close to 0. Both ^25j and ^20j consider only the 
case of -L(-) asymptotically constant and, for this case, their predictions differ. 

3.4. Relevance and irrelevance of the disorder: the results. The heuristic picture 
outlined in the previous subsection has now been made rigorous. A summary of these 
rigorous results is given in the next three theorems. We recall that /i^""(/3) = hc{0) — /9^/2 
and that the annealed free energy is f(0, h + /3^/2), so that the annealed critical behavior 
is obtained by looking at f(0, /i^""(/3) - (/3V2) + S) = f(0, hc{0) +6) as6\0. 

Theorem 3.5. Choose a £ (0, 1/2) and K(-) satisfying (jl.Sp and (jl.9p . Then there exists 
/3o > such that hc{(3) = h'^^^{(3) for (i < Pq. Moreover, for the same values of 13 the 
critical behavior of the quenched free energy coincides with the critical behavior of the 
annealed free energy: 

log f(/3, h,{P) + 6) ''^^ log F(0, hc{0) + 6). (3.25) 

Theorem [33] has been first proven in [4] by using a modified second moment method that 
we are going to outline in Section [3] It has been then proven also in [51], by interpolation 
techniques. Both works contain more detailed results than just Theorem 13.51 in particular 
(|3.25p has been established by showing that the stronger statement (j4.ip holds. As a 
matter of fact in [36] Th. 2.3] it has been proven that 

f(0,5)-f(/3,/i;^""(/3) + ,5) 



lim lim sup 

/3\0 5\0 



{(3y2){dsF{o,6)y 



0, (3.26) 
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written for hc{0) = for sake of compactness. Note that (|3.26|) is in agreement with (|3.17|) 
and in fact the first step in justifying that expansion. 

Theorem 3.6. Choose K{-) satisfying (jl.Sp and (|1.9p . If a > 1/2 we have h^P) > 
/i"""'(/3). Moreover 

(1) if a £ (1/2, 1) we have that for every e > there exists c^ > such that 

hc{0) - hT'^{(3) > c,/3^+^ (3.27) 

for every /3 < 1; 

(2) for every K{-) such if a > 1 there exists c > such that 

h,{P)-hT''{P) > cf3\ (3.28) 

for every (3 <1. 

The results in Theorem 13.61 are almost sharp, because in [H [5l] it is proven that for 
every K{-) such that a £ (1/2, 1) there exists C > such that 

/ic(/3) - /ir"(/?) < Cp^^, (3.29) 

for every /? < 1. On the other hand the bound in Theorem 13.6( 2) is already optimal (in 
the same sense) in view of the bounds summed up in the caption of Figure HI The result 
(j3.27p has now been improved to match precisely (j3.29p . i.e. it has been shown in [7] that 
in (j3.27p one can take e = and cq is still positive. 

Theorem 13.61 has been proven in [19] and we give an outline of the proof in Section [5l 
The method is based on estimating fractional moments of the free energy, while ()3.29p is 
derived by adapting the techniques yielding Theorem 13.51 (and a sketch of the proof is in 
Section[3|). In [19] the case q = 1 is not considered, but in [10] it is shown that the fractional 
moment method can be generalized to establish in particular that /ic(/3) > h'^^{P) also 
for Q = 1. 

For what concerns the critical behavior we have the following: 
Theorem 3.7. For every K{-) we have 



1 + Q 



\2 



(0 < F(/3, h) =) F(/3, h) - F{P, h,m < -^ {h - hc{P)r , (3.30) 

for every h (of course the result is non trivial only for h > hc{P)). 

The result in Theorem 13.71 has been established in [M] for rather general charge distri- 
bution. The proof that we give in Section [6] uses rather heavily the Gaussian character of 
the charges and it is close to the argument sketched in |33j . 

Let us point out that Theorem 13.71 coupled with Theorem 12.21 shows that the critical 
behavior of quenched and annealed systems differ as soon as a > 1/2, in full agreement 
with the Harris criterion: 

. ^ log (F(/3, h,{p) +5)- F(/?, K{P)))) a>l/2 log (F(0, K{0) + 5)- f(0, /l,(0)))) 

limmt > , 

(5\o log log 

(3.31) 

since the left-hand side is bounded below by 2, by Theorem 13.71 and the right-hand side 

is equal to max(l, 1/a), by Theorem 12.21 
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Remark 3.8. Theorem 13.71 therefore shows that the disorder, for a > 1/2, has a smooth- 
ing effect on the transition. The Harris criterion in principle is just suggesting that there 
is no reason to beheve that the critical behavior is the same. There is a general belief 
that disorder smooths the transitions: this is definitely the case for a number of statistical 
mechanics models to which a celebrated result of M. Aizenman and J. Wehr applies [3] 
(see also |38j). It should be however remarked that the Aizenman- Wehr smoothing mech- 
anism does not yield smoothing for the pinning model and that the argument leading to 
Theorem 13.71 is very different from the argument in [3j (for more on this issue see the 
caption of Figured]). 

Remark 3.9. The amount of smoothing proven by Theorem 13.71 was not fully expected. 
In fact in [48] it is claimed that for a > 1 the transition is still of first order, in disagreement 
for example with [17] [18] . Needless to say that it would be very interesting to understand 
what is really the value of the exponent for a > 1/2 and how it depends on a. In [5j it 
is shown that pinning models based on exponentially decaying inter-arrival laws may not 
exhibit smoothing. 

Remark 3.10. The results we have presented do not consider the case a = 1/2. The 
results in this case are incomplete and not conclusive under hypothesis (|1.9p . In [H[53] it is 
shown that under (|1.9p one has hc{(i) — h^^'^ {(5) < cexp(— l/(c/3^)) for some c > and j3 < 
1. This bound matches the prediction in [20]. It is not known whether hc{l3) — h'^^{(3) > 0, 
leaving open the possibility for the prediction in [25], i.e. hc{(3) = /i"""'(/3) for /3 small, to be 
the right one. One has to point out however that the approach in [25] is a development in 
powers of /3 that cannot capture contributions beyond all orders. This issue remains open 
and debated also at a heuristic level. About the critical behavior, the smoothing result 
in Theorem 13.71 is once again not conclusive (but it does imply for example that disorder 
smooths the transition as soon as lim„_^oo L(n) = 0, li L{-) is chosen as in Remark II. 4p . 
So, under assumption (II. 9p (which is the one that arises in the basic example, cf. ()1.4p ). 
the issue of whether a = 1/2 is marginally relevant or marginally irrelevant is open. 

4. Free energy lower bounds and irrelevant disorder estimates 

This section is mostly devoted to giving the main ideas of the proof of Theorem 13.51 
but in § 14.21 we will also explain why such arguments yield also (13.29P . 

4.1. The case of a < 1/2: the irrelevant disorder regime. As we have already 
pointed out, the annealed bound already yields, in full generality, that f(/5, h) < f(0, h + 
f3^/2), and hence hcif3) > hc{0) - 13'^ / 2 = h'^''''{(3). In order to pin, for a < 1/2, that 
hci(3) < /i""'"(/3) we need to prove a lower bound on the free energy showing that f(/3, h) > 
whenever f(0, h + /3^/2) > 0. We are actually aiming at capturing also the critical 
behavior of the free energy; in fact we are aiming at showing that for every e > there 
exists /3e > such that for every /3 G (0, /3e) we have 

liminf , ^'^2, , > 1-g, (4.1) 

h\o f(0, /i + /32/2) - ' ^ ^ 

which yields hc{f3) < h'^^{(3) and it is a result (sensibly) stronger than (j3.25p . Note that 
by what we have seen on the expansion of the free energy ()3.17p , or by the rigorous result 
()3.26p , we cannot aim at proving that the quenched free energy coincides with the annealed 
one. This actually casts some doubts about the applicability of second moment methods. 
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As a matter of fact if we choose h > — /3^/2 (/ic(0) =0), as usual with 5 = h + /3^/2 we 



can write 



varp ( Z 



'N,i 



^zi, 



^5 



exp 



I P 2^ InGrnr' J 



(4.2) 



Note the analogy with (|3.2U|) formula, in which 6 = 0: this time, since the underlying 
measure is P/i (see the beginning of Section[2]or Theorem 12 ■2p . rPlr' is a positive persistent 
renewal and the expression in (j4.2p is growing exponentially as A^ — > oo for every /? > 
(unlike for the 6 = case, in which the exponential growth sets up only for /? larger than 
a positive constant Pq). As we have pointed out, this was to be expected: can one still 
extract from ()4.2p some interesting information? The answer is positive, as shown by K. 
Alexander in f4j. 

The crucial point is not to take the limit in A^, but rather exploit (14. 2p up to the scale 
of the correlation length of the annealed system, which is just a homogeneous system with 
pinning potential 6 (see Remark 12. 4p . The idea is to establish that the quenched partition 
function is close to the annealed one with large probability up to the correlation length 
scale. Note that, as pointed out in Remark 12.41 o^i such a scale the annealed partition 
function starts exhibiting exponential growth and one feels the size of the free energy (we 
will actually need to choose the size of the system, call it Nq, to be a large, but finite, 
multiple of the correlation length, the relative error e in (j4.ip leaves the room to make 
such an estimate). Once such an estimate on a system of length A^o is achieved, it is a 
matter of chopping the polymer into N/Nq portions: some work at the boundary of these 
regions is needed and for that we refer to [1] , while we focus on explaining why the second 
moment method works up to the correlation length scale. 

Let us therefore go back to (j4.2p and let us set Nq := q/F{0, 6) (and assume that it is an 
integer number). With such a choice, by exploiting the estimates outlined in Remark [ 
it is not difficult to see that for every K{-) and q > there exists cxiq) > such that 



N 



K^^o,- 



Eexp 6^1neT < CK{q)- 



(4.3) 



n=l 



This is because, as long as q is finite and 6 tends to zero, the system is in the critical 
window (as a matter of fact, in [50] the limit of EZj^ ui ^^ ^ tends to zero, q kept fixed, is 
computed and the notion of critical window is further elaborated). The constant cxiq) of 
course diverges as g /" oo. On the other hand EZ^^ > 1 so that in order to estimate the 
quantity in (j4.2p it suffices to estimate 



E^ 



N 



N 



exp 



(5^(l„e^ + l 



n€T' J 



exp 



P^^l 



nGrHr' 



1 



n=l 



n=\ 



(4.4) 



We now use the Cauchy-Schwarz inequality and the fact that (exp(x) — 1)^ < exp(22;) — 1 
for x > to bound the expression in ()4.4p (and therefore the expression in ()4.2p ) by 



E 



expf2(5^1„eJ E ( exp ( 2/3^ ^ l„e,nr' j - 1 j 



=: Ti • Ta. 



(4.5) 



But, by ()4.3p . T\ is bounded by a constant (which depends on q). On the other hand, T2 
has been already estimated in (I3.20p - (l3.24p . and it is 0(/?), thanks to the fact that the 
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renewal r n r' is terminating (since a < 1/2). Therefore, by Chebychev inequality, for 
every e > 

P (Z^o,. > (1 - 6)EZ^„, J < ^^V?/?, (4.6) 



where Ck{q) is a constant depending on K{-) and q (it is just the constant cxiq) of 
when 5 is replaced by 26) and c is taken form (|3.24p . Since EZ^ ^ is bounded below by 
exp(— (5) exp(F(0, 6)No) = exp(— (5 + q) (see ()2.18p ) we see that on the scale of correlation 
length the quenched partition function grows (almost) like the annealed one with large 
probability if /? is small enough. 

4.2. Lower bounds on the free energy beyond the irrelevant disorder regime. 

The technique for lower bounds on the free energy that we have outlined, as well as the 
technique in [Ml [56]) lead to upper bounds on hdP) also in the case a G [1/2, 1), see ()3.29p 
and Remark 13.101 But a look at Theorem I3.6l f2l suffices to see that the case q > 1 is 
somewhat different, because it is no longer true that hc{l3) — h'^"-^{j3) = o{j3'^). So, in this 
case, the easy bound /ic(/3) < ^c(O) (which is just a consequence of convexity) is optimal 
in the sense that hc{P) - /i^""(/3) < 13"^ / 2 = 0{p^). 

Let us therefore explain why the second moment method yields also (j3.29p when a G 
(1/2, 1). For this we go back to ()4.2p and (j4.4p (we are still placing ourselves on the scale of 
the correlation length). The term Ti is still bounded by a constant, that of course depends 
on K{-) and q, just as in the a < 1/2 case. The term T2 this time grows exponentially in 
A^, because this time r n r' is persistent, and we have to worry about the size of A^ also 
for this term. But let us quickly estimate the growth rate of T2 and for which values of A^ 
we can expect this term to be small for /? small. A necessary condition, that with some 
careful work one can show also to be sufficient, is that the expectation of the term in the 
exponent in the exponent is small, namely that (cf. Proposition II. 2p 

N 
/32^P(n e rf ^^"^ c,/3'A^'°-\ (4.7) 

n=l 

has to be chosen small. However we still keep N = Nq = q/F{0, 6), that is A^ of the order 
of (5~^'" (by Theorem 12. 2p . times a constant which is large if q is large. Plugging such a 
value of A^ in ()4.7p we see that we are asking p'^S^^~'^°'''^'^"' to be small. Therefore in this 
regime we expect the second moment method to work, leading to localization and also to 
the fact that the quenched free energy is fairly close to the annealed one, if 

6 > c/32°/(2°-i), (4.8) 

with c a small (fixed) constant. But this what is claimed in (I3.29p . 

5. Relevant disorder estimates: critical point shift 

Annealing is the standard procedure to get upper bounds on disordered partition func- 
tions. One can go beyond by partial annealing procedures, like the constrained annealing 
procedure [44J, and this does give some results, see e.g. P], but for pinning models con- 
strained annealing, in the infinite volume limit, yields nothing beyond the annealed bound 
if we are concerned with identifying the critical point [15]. There is therefore the need for 
a different idea. 
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5.1. Fractional moment estimates. A tool that allows to go beyond the annealed 
bound hciP) > h'l^^{j3) (we set /ic(0) = also in this section) turns out to be estimating 
Ai^ := E[(Z^^)'''] for 7 G (0, 1) by means of the basic inequality 

7 



Y^aA <Y.a], (5.1) 



that holds whenever Uj > for every j. This has been pointed out by F. L. Toninelli in 
[55| . Inequality (j5.ip has been exploited also in other contexts, notably in |221 I14j. to get 
upper bounds on the partition function of the directed polymer in random environment, 
and in [^ to establish localization of eigenfunctions for random operators, in particular in 
the Anderson localization context. 

For pinning models one applies ()5.ip to the renewal identity 

N~l 

Zn,u. = Y. ^n,u.K{N - n)e7v, with Cn := exp(/3c^^ + h), (5.2) 

n=0 

and, by taking the expectation, one gets to 

N-l N 

An <^ [Ci] Y ^n^(^ -^y = Y1 ^N-nQ{n), (5.3) 

n=0 n=l 

where Qn := K[S,i]K (n)'^ . Now the point is that (j5.3p implies 

An < (X^Ql?^) I _ max _ An, (5.4) 

\n=l J "-O'l'-'^ 1 

SO that, if ^^ Qin) < 1 we have An ^ Aq = 1 for every A^. Summing everything up 

oo 

HCi]y"K{ny < 1 =^ supAn<1. (5.5) 

And of course if An has sub-exponential growth the free energy is zero since 

^ElogZ^,^ = -LElog(Z^J" < ^logAN. (5.6) 

Remark 5.1. The discrete convolution inequality (|5.5p can actually be exploited more. 
Observe in fact that the solution to the renewal equation 

Af 

uq = 1 and UN = /_] Q('^)^Ar-n for A^ = 1,2,..., (5.7) 

n=l 

dominates An- But if Q{-) is a probability distribution (possibly adding Q{oo)), then u. 
is the renewal function of the (5(-)-renewal and if YlnQi^) < -*- then un ~ cQ{N) for 
N —* oo (c is an explicit constant, see Theorem 1 1.2p . Therefore 

oo 

E [^7] Y K{ny < 1 =^ there exists C>0 such that An < CK{Ny. (5.8) 

n=l 
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We are now left with verifying for which values of (3 and h we can find 7 such that 
E[^^] Yl'^=i K{n)'^ < 1. Let us consider the case ^„ K{n) = 1: in this case Yin K{n)'^ > 1, 
so the question is for which values of (3 and h the pre-factor E^^ is sufficiently small. This 
is a straightforward computation [5 = h + /3^/2): 

EC? = exp (- y7(1 - 7) + 'J7) , (5.9) 

which is small for (5 sufficiently large, for every fixed value of 6. This result is therefore 
saying that (it may be helpful to keep in mind Figure H]): 

(1) K{I3) > h'^'^'^ip) if P is sufficiently large; 

(2) the gap between hc{P) and /i""'"(/?) becomes arbitrarily large as /? tends to infinity: 
in fact it is of the order of /3^ . 

Note that this approach yields very explicit bounds, but it does not give results for small 
values of (3. 

5.2. Iterated fractional moment estimates. To go beyond the estimate we have just 
presented, in [19j another renewal identity has been exploited, namely: for every fixed k 
and every N > k 

N k-1 

Zh,u. = E ^N-n,. E ^(^ - ^■) ^N-jZlg^^,^. (5.10) 

n=k j=0 

This is simply obtained by decomposing the constrained partition function according to 
the value N—n of the last point of r before or at N—k (0 < N—n < N—k in the sum), and 
to the value N — j of the first point of r to the right of N — k (so that N — k < N — j < N) . 
Of course Z^^jy-j^ has the same law as Z^^ and the three random variables Z^_^ ^, S,n~j 
and Z'^gj^^j^ are independent, if n > fe and j < k. 

N-n N~k N -j 



X H — X X X-H X X-H H X • X — X X 

^ ZN-n,u, K{n-j)^N-j ^Af-i,Af,c^ ^ 

Figure 5. The renewal identity (|5.10p is obtained by fixing a value of k and summing 
over the values of the last contact before N — k (the large dot in the figure) and the first 
contact after N — k. The two contacts are respectively N — n and N — j (crosses are 
contacts in the figure). 

Let < 7 < 1, set once again An := E[{Zf^ ^y] and use (|5.ip and (|5.10p to get for 
N>k 

N k-1 

An < nCi] E ^^- E ^(^ - ^y^r (5-11) 

n=k j=0 

This is still a renewal type inequality since it can be rewritten as 

N 

An < E^iV-nQfc(n), (5.12) 

n=l 
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with Qk{n) := IE[^J'] 'Y1,j=q K{n — j)^Aj if n > A: and Qk{n) := for n < k. In particular 
if for given (3 and h one can find A; S N and 7 G (0,1) such that 

00 k—l 

p--=Yl Qk{n) = nCi] E E ^(^ - ^r^j < ^ (5-13) 

n n=fc j=0 

then one directly extracts from (j5.12p that 

An < praax{Ao,...,A]y_k}, (5.14) 

for N > k, which implies that An < max{j4o, . . . , ^fc-i} and hence f(/?, /i) = 0. 

Remark 5.2. Like in Remark 15. II one can be sharper by exploiting the renewal structure 
in (j5.12p . The difference with Remark 15.11 is that in this case N > k. In order to put 
()5.12p into a more customary renewal form we set An := A]y'i-N>k, so that 

N-{k-l) ^ fc-1 

An < Yl AN-nQkin)+PkiN), with PkiN)=Y,Mn)QkiN-n), (5.15) 

71=1 n=0 

and therefore there exists c > (depending on k, K{-) and 7, besides of course /? and h) 
such that Pk{N) < cQk{N). Let us now consider the standard renewal equation for the 
(5fc(n)-renewal: uq = 1 (but of course one can choose an arbitrary uq > 0) and 

N N-{k-l) 

UN = '^UN-nQkin) = E UN-nQk{n)+UoQk{N), (5.16) 

n=l n=l 

where the first equality holds for A^ = 1,2,... and for the second one we have used that, 
since Qk{n) = up to n = A; — 1, we have ui = U2 = ■ ■ ■ = u^-i = 0. Once again if 
Ylin Qkin) < 1, that is if /9 < 1, we are dealing with a renewal equation of a terminating 
process and therefore un behaves asymptotically like (a constant times) Qk{N). Compar- 
ing (j5.15p and (j5.16p one obtains that there exists a constant C = C{K{-),k,^,h, (3) such 
that 

An < CK{Ny. (5.17) 

5.3. Finite size estimates by shifting. What the iterated fractional moment has done 
for us is reducing the problem of estimating the free energy from above to a finite volume 
estimate. Notice in fact that estimating p, cf. ()5.13p . amounts to estimating only (a 
fractional moment of) Z^^, for j < k (note the parallel with Remark 13.41 ). This type 
of estimates demands a new ingredient, which is more easily explained when a > 1. A 
preliminary observation is that p is bounded above by e~^ 'l2i=oi^j/(^ — j)(i+")7-i^ with 
e a constant that depends on K{-) and 7 so that p < 1 is implied if for a given 7 

^ A- 

Note that, unlike the case treated in § 15. H here the pre- factor E^^, and therefore /3 and 
h, has only a marginal role: as long as /? and h are chosen in a compact, which is what 
we are doing since we are focusing on the critical region of the annealed model at small 
or moderate values of /3, e can be chosen independent of the value of /3 and h. We shall 
see that the expression in (I5.18P can be made small for example by choosing k large. 
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Let us start by observing that we know of course (Jensen inequality) that for h 



^0 < (E^,.)' 



E 



exp I 5 ^ Iner I ; i G r 



exp(7F(0,(^)j)P5(j G r)" 



(5.19) 

We are of course interested in 5 > 0: by the Renewal Theorem Y*s{j £ t) is bounded below 
by a positive constant (even if 5 were zero!), so this term cannot be of much help and we 
simply bound it above by one. On the other hand the exponentially growing term stays 
bounded for j up to the correlation length of the annealed system {cf. Remark 12. 4p : we 
therefore choose k := 1/f(0, (J) (again, assume that it is in N). At this point we observe 
that we can choose 7 G (0, 1) such that 



(l + a)7 > 2, 



(5.20) 



then the expression in (j5.18p is bounded for k large, that is 5 small. This is not yet what 
we want, but a more attentive analysis shows that one has 



k-l-R 

E 

j=0 



A 



(fc_j)(l+a)7 



^^:^^y^<exp(7)5:r(^+")^+^ <e/2, 
j>R 



(5.21) 



for any k > R and R chosen sufficiently large (depending only on 7, a and e). This has 
been achieved by using (j5.19p . We have therefore to show that 



fc-i 



Ai 



^.^^^(,_,-)(i^.).- 



T ^ e/2- 



For this we set 



If we are able to show that 



A, 



lim sup max Aj . 

5\0 j=k-R,...,k~l 



R 



^^E 



-((l+a)7-l) 



< e/3, 



(5.22) 
(5.23) 

(5.24) 



i=l 



then (j5.18p would be established (for 5 small and k = 1/f(0, 5). Of course in (j5.24p one can 
replace R with 00 obtaining thus a more stringent condition (but, in the end, equivalent, 
since we are not tracking the consants). For a proof of (j5.24p one has to go beyond (|5.19p 
and in doing so the size of /3 turns out to play a role. 

In order to go beyond (|5.19p the new idea is a tilting procedure (first proposed in [31j), 
that, given the Gaussian context, reduces to a shift. The idea is based on the following 
consequence of Holder inequality 

1-7 



A, 



E' 



(^J 



dP' 



M 



< 



E' [Z. 



17 TIT./ 



J,^J 



E' 



dP 
dF 



[to 



1/(1-7)' 



(5.25) 



where P' is a probability with respect to which P is absolutely continuous. In order to 
make the choice of P' let us fix (5 = a/3^, o a constant that we are going to choose along 
the way: P' is the law of the sequence 



LOl 



V^, 



OJ2 



/a/32 



,Wfc — '\/a(3'^,u;k+i,'^k+2, ■ ■ ■ 



(5.26) 
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which is a sequence of independent (non identicahy distributed) variables, 
readily computes 



E' 



_dl 
df 



[UJ 



l/(l-7)' 



1-7 



exp 



7 



1-7 



al3^k 



One then 



(5.27) 



but a0^k = 5/f(0, 5) and the ratio (5/f(0, 5) tends to a positive constant as 5 \ since 
a > 1, cf. Theorem 12.21 Let us now turn our attention to E'Z? which, for j < k, 
coincides with K'Z'^ , — -. But this is just the partition function of a homogeneous 

model with negative pinning potential if we choose a small, namely (for conciseness we 
look only at the case j = k) 



^'ZL 



E 



E 



exp 



exp 



I -/3^(\/a - a) X] 1«G^ \;k£T 

)i|:w):*€ 



6 \ 1 



F(0,<5), 



(5.28) 



But limfc_^oo(l/^) X]n=i l^er = 1/E[ri] P-a.s. and this readily implies that WJZ^^ is 
bounded by a constant that can be chosen arbitrarily small, provided one chooses a suffi- 
ciently small (so a~^'^ — 1 is large). The argument easily extends to j between k — R and 
k so that (j5.24p is proven since A]^ can be chosen arbitrarily small for a sufficiently small 
and every 5 < 5q (for some 5q > 0). This concludes the argument for the case a > 1, that 
is Theorem 13. 6r 2). 

The case a G (1/2, 1) (Theorem 13.6( 1)) is conceptually not very different: the main 
difference lies in the fact that it is no longer sufficient to show that Aj is small for j close 
to k, one has actually to extract some decay in j. But the fact that Aj does decay with 
J, at least if j < 1/f(0,5), is already rather evident from (J5.19P from the fact that the 
term Y's{j £ r) is, at least till j < k = 1/f(0, (5), close to P(j G r) which behaves for j 
large as j~(^~") (times a constant, cf. Theorem ll.2p . The argument is however somewhat 
technical and we refer to [I9j for details. 



6. Relevant disorder estimates: the critical exponent 

The argument leading to Theorem 13.71 is based on the rare stretch strategy sketched in 
Figure [6l It is based on a one-step coarse graining of the environment on the scale £ G N, 
I ^ ^ ^ A^. Actually one should think of i as very large but finite. We assume N/i G N 
and we look at the sequence of IID random variables defined as 



Yj 



1e^, with Ej = jw : log ^^"^go-iK^ > aF{(3, /i + 5)} , 



(6.1) 



where 5 > 0, a G (0, 1) (eventually a y 1) and j = 1,2,..., but of course only the j's 
up to N/i are relevant to us. Note that the Y variables are Bernoulli random variables 
of parameter p{i) := P(i?i) and, since a < 1, p{l) is small when i is large, by the very 
definition of the free energy and its self-averaging property {cf. Theorem 13. ip . One can 
actually show rather easily that 

'2/52' 



liminf — \og p{i) 



> 



(6.2) 
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We give a proof of this inequality below, but the intuitive reason is that the probability 
of observing J2i=i^i ~ ^S/P behaves like exp(— M^/2/?^) for £ large (this is the standard 
Cramer Large Deviation result). When such a Large Deviation event occurs, the environ- 
ment in the £-block will look like the original to variables translated of S/P, that is /3uji + h 
looks like /Jwj + h + S. And in that block the logarithm of the partition function will hence 
be close to ^f(/3, h + 6). Shifting the mean is of course only one possible strategy to make 
the event Ei typical and hence such an argument yields only a lower bound on p{i). 




^3^ " gJ' 



Figure 6. The rare stretch strategy is implemented by looking at blocks of ui variables 
of size £. To block j is associated a Bernoulli random variable Yj : such a random variable 
is a function of the ui variables in the block and it determines whether a system of the size 
of that block, with constrained boundary conditions and precisely with the ui variables 
of the block, has a sufficiently large (in fact, an atypically large) partition function. The 
precise definition of Yj is in (|6.1|l . In the figure Y4 = Yn = Y12 = 1, while the other 
Y variables are zero. The gaps between the success blocks are parametrized as Gn^, so 
that {G„}„ are IID Geometric random variables: P(Gi — k) = [1 — p)''p (fc = 0, 1, 2, . . . 
and p = p{i), given in the text). The lower bound is then achieved by restricting the 
partition function to the r trajectories that visit only success blocks and visit the first 
and the last point of the block (also if two success blocks are contiguous: this is the 
case of the second and third success blocks in the figure) . Such a strategy is profoundly 
different from that employed in [3] , that is based on the effect of typical fluctuations (on 
the Central Limit Theorem scale) and competition with boundary effects. Our strategy 
is instead a Large Deviation strategy and, in a sense, it exploits the flexibility of the 
polymer to target rare regions (the boundary conditions play no role). 



We now make a lower bound on Zj^^ by considering only the r-trajectories that visit 
all and only the i blocks for which Yj{uj) = 1 (see Figure [6]). The renewal property leads 
to a rather explicit lower bound, namely (with the notation of the figure) 

logZ^,^ > 2^ log^.Vi)., + X^ logi^(G,£) + 0(logiV), (6.3) 



j<N/l: 



n=\ 



where the 0(log N) term comes from the last excursion. Let us now divide by N and 
take the limit, keeping into account that, by definition of the Y variables, we have a lower 
bound on the partition functions logZ^ (^j_^i that appear in the right-hand side. We 
therefore obtain 



A^^oo N N^oo \ N J\Iy{uJ) 



Y, log K{GJ 



n=l 



p(i) 

ap{£)F{f3, h + 6) + ^-V" limsup 



Afvii^) 



i 



N^CXD 



Mriuj) 



Y. logK{Gr. 



n=l 



(6.4) 
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where in the last step we have used the strong law of large numbers (the limits are in 
the P(dti;)-a.s. sense) to estimate the leading behavior of the number of successes in 
an array of N/(. Bernoulli variables of parameter p{t). Moreover, to be precise, (|6.4|) 
holds if we set K{G) = 1 (which we do only here). The (superior) limit that is left in 
the expression is also easily evaluated by using the strong law of large numbers after 
having observed that, since \ogK{x) ~ —(1 + a)logx, when G„ = 1,2,... we have 
\og K{Gnt) > —a~^{l + a)(logG„ + log£) for £ sufficiently large (uniformly in the value 
of Gn'- note that 1/a is once again just a number larger than 1). The net outcome is 
therefore 

F(/?,/i) > ap(^)F(/3,/i + 5)-^a-Hl + a)(E[logGi;Gi >0]+P(Gi >0)log^), (6.5) 

Since Gi is a geometric variable of parameter p{i) we directly compute E[logGi;Gi > 
0] = (1 + o^(l))log(l/p(^)) which, by dO]), is bounded above by a'^ 15"^ / {2(3'^) (i large: 
once again a^^ is just used as an arbitrary constant larger than one). Therefore 



FiP,h) >p{l) 



2/32 



aF(/3, h + 5)- a-^ \J.^ + 0^(1) 



(6.6) 



where the term 0^(1) is —c{\ogt}/l (c > 0) and this bound holds, given a G (0,1), for 
every i larger than some Iq. 

Now we set h = hdfi) in (j6.6p . so that the left-hand side is zero and therefore 

2/32 

for every I > Iq, so that f(/3,/ic(/3) + 6) < a'^ {6'^ {I + a)/ (2/3^)), and since a £ (0,1) is 
arbitrary we can let a ^ 1 and we are done. 



aF(/3, h,iP) + 5)- a-^ ^ , ' + o,(l) < 0, (6.7) 



For completeness we give a proof of (|6.2|) . We call P^ the law of the sequence of random 
variables 

uji +5/f3,L02 + 6/13, . . . ,LOi + 5/l3,u;i+i,iOi+2, ■ ■ ■ (6.8) 

Note that ¥i{Ei) tends to one as i becomes large, simply by definition of free energy and 
because the law of {(3uJn-\-h}n=i,.../-, when uj is distributed according to P^, coincides with 
the law of {f3ijJn + h + 5}n=i,...,i, when u is distributed according to P. We compute the 
relative entropy 

-, r Aii. 1 .ir2 

(6.9) 



^ + - , (6.10) 



n (v^ p") := E^ 


"1 dP£, " 

log ^p {^) 


2P^ ' 


and by a standard entropy inequality (see e.g. [29, § A. 2]) 


log .^^^ > -^ W(P£P)+- = 


1 (5 


FfXEi) W 


and this vields (j6.2|). 
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